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Abstract 

A well-known diffuse interface model consists of the Navier-Stokes equations 
nonlinearly coupled with a convective Cahn-Hilliard type equation. This system 
describes the evolution of an incompressible isothermal mixture of binary fluids and 
it has been investigated by many authors. Here we consider a variant of this model 
where the standard Cahn-Hilliard equation is replaced by its nonlocal version. More 
precisely, the gradient term in the free energy functional is replaced by a spatial con- 
volution operator acting on the order parameter tp, while the potential F may have 
any polynomial growth. Therefore the coupling with the Navier-Stokes equations is 
difficult to handle even in two spatial dimensions because of the lack of regularity 
of ip. We establish the global existence of a weak solution. In the two-dimensional 
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case we also prove that such a solution satisfies the energy identity and a dissipative 
estimate, provided that F fulfills a suitable coercivity condition. 



Keywords: Navier-Stokes equations, nonlocal Cahn-Hilliard equations, incom- 
pressible binary fluids, existence of weak solutions. 
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1 Introduction 

A well-known model which describes the evolution of an incompressible isothermal mixture 
of two immiscible fluids is the so-called model H (see [27J [25], cf. also [HI [301 E2] and 
references therein). This is a diffuse- interface model (cf. [I]) in which the sharp interface 
separating the two fluids (e.g., oil and water) is replaced by a diffuse one by introducing 
an order parameter ip. The dynamics of tp, which represents the (relative) concentration 
of one of the fluids (or the difference of the two concentrations), is governed by a Cahn- 
Hilliard type equation with a transport term. This parameter influences the (average) 
fluid velocity u through a capillarity force (called Korteweg force) proportional to /iVy?, 
where \x is the chemical potential. Note that this force is concentrated close to the diffuse 
interface. 

In a simplified setting where the density g of the mixture is supposed to be constant 
as well as the viscosity v and the mobility m, the model reduces to 



in Q x (0,T), where Q is a domain in M d , d = 2, 3, T > is a given final time, tt is 
the pressure, k is a given positive constant and h represents volume forces applied to the 
binary mixture fluid. The chemical potential fi is the first variation of the free energy 
functional (see [IB]) 



Here F represents the (density of) potential energy. This function is usually a double- 
well potential whose wells are located in the pure phases, while £ and r\ are given positive 
constants. The potential can be defined either on the whole real line (smooth potential) 
or on a bounded interval (singular potential). The latter case (in a logarithmic form) is 
the most appropriate choice from the modeling viewpoint (cf. [IB]), while the former can 
be considered as an approximation. 



ip t + u ■ V<p = mA/i 

QU t — uAu + (u ■ V)li + V7T = KflVf + h 

div(n) = 



(1.1) 
(1.2) 
(1.3) 




(1.4) 
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In the context of statistical mechanics, the square gradient term in (11. 4p arises from 
attractive long-ranged interactions between the molecules of the fluid and £ can be related 
to the pair correlation function (see, e.g., [4] and references therein). We also recall that 
K and £ are of the same order as the interface thickness e > 0, while rj is proportional to 
e _1 . On account of (ll.4p . the chemical potential takes the following form 

fi = -ZAp + V F'(tp). (1.5) 

Systems like (ll.ip ~ fll.5p . also known as Cahn-Hilliard-Navier-Stokes systems, have been 
studied from the mathematical viewpoint by several authors (see, for instance, [TJ El |3j 
EH ED H21 ED E21 E2 EH], cf. also H3 dSl US] for numerical issues). 

A different form of the free energy has been proposed in [231 EI] and rigorously justified 
as a macroscopic limit of microscopic phase segregation models with particle conserving 
dynamics (see also [2]). In this case the gradient term is replaced by a nonlocal spatial 
operator, namely, 



= 7 / / J(x-y)((p(x)-(p(y)) 2 dxdy + Tj I F((p(x))dx, (1.6) 
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where J : M d — > K is a smooth function such that J(x) = J(—x). Taking the first variation 
of £ we can define the chemical potential associated with the nonlocal model 

fi = cup — J * (p + r}F'(ip) (1.7) 

where 



(J * ifi)(x) :— / J(x-y)(p(y)dy, a(x) := / J(x-y)dy, x G fi. (1.8) 

The corresponding nonlocal Cahn-Hilliard equation (pt = mA/j, can be derived from ideal- 
ized microscopic models through suitable limits like the diffusion equation and the Boltz- 
mann equation. Moreover, the evolution in the sharp interface limits are the same as those 
derived from the classical Cahn-Hilliard equation in the corresponding limits (see [24J ) . 
However, from the mathematical viewpoint, the nonlocal Cahn-Hilliard equation, due to 
its integrodifferential nature, is rather difficult to handle (see, e.g., [T|. 18], [15 | [T9 |. [20 ]. [26 ], 129] ) . 
Here we consider system (ll.ip -( fO|) with (JTTTJ) . More precisely, taking for simplicity all 
the constants but v equal to one, we want to study the following initial and boundary 
value problem 

ipt + u-Vp = A/2 (1.9) 
fi = a<p- J*<p + F'(<p) (1.10) 
u t - vAu + [u ■ V)it + Vvr = /iVy? + h (1.11) 
div(w) = (1.12) 



3 



-C = 0, u = on<9fix(0,T) 
on 

u(0) = Uq, <p(0) = ipo in Q, 



(1.13) 
(1.14) 



where Q C M. d , d = 2, 3, is a bounded domain with sufficiently smooth boundary and 
unit outward normal n. The no-flux boundary condition for \i is the usual one for Cahn- 
Hilliard type equations (cf., e.g., [7]) and implies the conservation of mass (see Remark 
[5] below). The no-slip boundary condition for u is also standard especially when one 
wants to investigate a new model involving Navier-Stokes equations (periodic boundary 
conditions can also be considered). 

In this contribution we prove the existence of a global weak solution for smooth po- 
tentials F of arbitrary polynomial growth. Moreover, if F satisfies a suitable coercivity 
condition then we can slightly improve the smoothness properties of the solution. In 
particular, we show the validity of an energy identity if d — 2. These results are a first 
step towards the mathematical analysis of problem fll.9l) - fll.l4l) . However, further issues 
(such as, e.g., uniqueness in two dimensions) do not seem so straightforward to prove. 
The main difficulty arises from the presence of the nonlocal term which implies that (p is 
not as regular as for the standard (local) Cahn-Hilliard-Navier-Stokes system (cf. Remark 
M below). For this reason, we have not been able even to establish uniqueness of weak 
solutions in two dimensions. 



2 Notation and functional setup 

Let us set V s := D{B S / 2 ) for every seK, where B = — A + 1 with homogeneous Neumann 
boundary conditions. Hence we have 

V 2 = D(B) = Ive H 2 (VL) : — = on dfl L 

We also define H := Vq = L 2 (Q) and V :— V\ — if 1 (fi). Then we introduce the classical 
Hilbert spaces for the Navier-Stokes equations (see, e.g., [35J) 

G dlv := {u e CS°(Q) d : div(u) = Of 2 ^, 

and 

V div := {u e Hl{n) d : div(u) = 0}. 

We denote by || • || and (•, •) the norm and the scalar product, respectively, on both H and 
Gdiv We recall that Vdiv is endowed with the scalar product 

(u, v) Vdiv = (Vu, Vt> ), Vtt, v E V div . 
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We also need to introduce the Stokes operator A : D(A) D Gdiv — > Gdiv Recall that, in 
the case of no-slip boundary condition f )1.13p 

A = -PA, D(A) = H 2 (n) d nv dlv , 

where P : L 2 (Q) d — >■ Gdiv is the Leray projector. Notice that we have 

{Au, v) = (u, v) Vdiv = (Vm, Vv), Vm g D(A), W G 

We also recall that A^ 1 : — > Gdiv is a self-adjoint compact operator in Gdiv and by 
the classical spectral theorems there exists a sequence Xj with 

< Ai < X 2 < ■ • ■ , Xj -»■ oo, 

and a family of u>j G -D(^4) which is ortho normal in Gdiv and such that 

Aifj = XjWj. 

Finally, for u, v, w G we define the trilinear V^,— continuous form 

b(u,v,w) = / (w • V)t> • u>, 
in 

and the bilinear operator B from Vdi V x into defined by 

(23(u, f ), w) = w), Wu, v, w G Vrfj„. 

We recall that we have 

= -b(u,v,w), Vu,v,weV d iv, (2.1) 

and that, for every u, v and u> G V^,, the following estimates hold 

\b(u,v,w)\<c\\u\\ 1/2 \\Vu\\ 1/2 \\Vv\\\\S/w\\, fore/ = 3, (2.2) 
|6(u,v,w)| < cllnll^llVnll^UVwIlll^ll^llV^II 1 / 2 , for d = 2. (2.3) 

In this paper c will stand for a nonnegative constant depending possibly only on J, 
f,Q,v and T. The value of c may vary even within the same line. We shall denote by 
N, M or L generic nonnegative constants that depend on the initial data uq, ipo and on 
h and whose values will be explicitly pointed out if needed. 

3 Main result 

In this section we first define the notion of weak solution to problem fll.9p - fll.14p which will 
be called Problem P. Then we state the main result of this paper and a related corollary. 

Our assumptions on the kernel J, the potential F and the forcing term h are the 
following (cf. also (11.81) ) 
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(HI) J G W 1,1 (R d ), J{x) = J{-x), a(x) := / J(x - y)dy > 0, a.e. igSI. 

Jn 

(H2) F G C 2 (IR) and there exists c > such that 

F"0) + a(x) > c , Vs G R, a.e. x6H. 

(H3) There exist ci > |||J||L 1 (K d ) an d c 2 G IR such that 

F(s) > ds 2 - c 2 , Vs G E. 

(H4) There exist c 3 > 0, c 4 > and p G (1, 2] such that 

\F'(s)\ p <c 3 \F(s)\ + c 4 , VsGM. 

(H5) / i GLL(0,T;yj i JforallT>0. 

Remark 1. The requirements of assumption (HI) are standard for the nonlocal Cahn- 
Hilliard equation (see, e.g., [7] for slightly stronger hypotheses). 

Remark 2. Assumption (H2) implies that the potential F is a quadratic perturbation of 
a (strictly) convex function. Indeed, if we set a* :— ||a||oo> then F can be represented as 

F(s) = G(s)-^s\ (3.1) 
with G G C 2 (1R) strictly convex, since G" > Cq in Q. 

Remark 3. Assumption (H4) is fulfilled by a potential of arbitrary polynomial growth. 
In particular, (H2)-(H4) are satisfied for the case of the physically relevant double-well 
potential, i.e. 

F( S ) = (1-, 2 ) 2 . 

In this case we take p = 4/3 in (H4), while assumption (H2) is satisfied if and only if we 
have p > Co + m , where m = — min se R F"(s). 

By weak solution we mean 

Definition 1. Let u G G<nv, y?o £ H with F(ip ) G L x (f2) and < T < +oo be given. 
Then [u, ip] is a weak solution to Problem P on [0, T] corresponding to uq and (po if 

• u, (p and fi satisfy 

u G L°°(0,T;G dlv )nL 2 (0,T;V dlv ), (3.2) 

u t eL^(0,T;V^ lv ), if d = 3, (3.3) 

M( GL^(0,T;l/ rf 'J, V 7 G(0,1), if d = 2, (3.4) 
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p G L°°(0, T; H) fl L 2 (0, T; V), (3.5) 

^eL 4/3 (0,T;V"), i/ d=3, (3.6) 

V9 t GL 2 - 5 (0,T ; y'), We (0,1), if d = 2, (3.7) 

pGL 2 (0,T;V). (3.8) 



setting 



p(x,<p) :=a{x)tp + F'{y), (3.9) 
then, for every ip e V, every v € and for almost any t e (0, T) we /iave 

(<^, V) + (Vp, V?A) = / (u ■ V*p)p + f{VJ*p)- W, (3.10) 

v) + i/(Vtt, Vv) + &(«, u, f ) = - / (v • V/i)v? + (/i, u). (3.11) 



• the following initial conditions hold 

u{0)=u , (p(0) = (p o . (3.12) 

Remark 4. Since p = p + J * p, from Definition [1] we have that p e L 2 (0,T; I 7 ). 

Remark 5. ft is immediate to see that the total mass is conserved. Indeed, choosing 
ip — 1 in f)3.10p . we have (<p t , 1) = whence (<p(i), 1) = (<po, 1) for all t > 0. 

Remark 6. The initial conditions f)3.12p are meant in the weak sense, i.e., for every v e 
we have («(£), i>) — >■ ("Uo> f ) as t — )• 0, and for every x e V we have (<£>(t), x) ~^ (^Qj x) 
as t — )• 0. ft can be proved that w € C TO ([0, T]; GWu) and p G C w ([0, T]; if). 

Theorem 1. Let u G Gdi V , <Po G if suc/i i/jcrf F(po) G L 1 (f2) and suppose that (H1)-(H5) 
are satisfied. Then, for every T > t/iere exists a u>ea£; solution [u, <p\ to Problem P on 
[0, T] corresponding to Uq, (po with p t satisfying 

V t EL°°(0,T-X), if 1<p<1 ± T , s= ^- d ^+ 2d , 

if, 6 L~(0,T; O n L'(0,T; i/ p = s > r > 2, 

2 a — 1 z 

^ G L 2 P /( 2P -3) (0;T . ^ z/ rf = 3; 3/2 < p < 2, S = 
Furthermore, setting 

£{u{t),p{t)) = \\\u{t)\\ 2 + - [ [ J(x-y)(p(x,t)-p(y,t)) 2 dxdy + [ F(<p(t)) 
A 4 Jn Jn Jn 

the following energy inequality holds for almost any t > 

S(u(t),p(t)) + [\u\\Vu(T)\\ 2 +\\Vp(r)\\ 2 )dr<£(u ,p )+ [ (h(r),u(r))dT. (3.13) 
Jo Jo 
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On account of the typical examples of double- well smooth potentials (cf. Remark |3]), 
the following additional assumption sounds reasonable (see, e.g., [7, (A2)]) 



(H6) F G C 2 (IR) and there exist c 5 > 0, c 6 > and q > such that 

F"(s) + a(x) > c 5 \s\ 2q - c 6 , VsGM, a.e. x G O. 

This requirement can replace (H3) in the proof of Theorem [T] (see (13.141) below). Indeed, 
(H6) implies the existence of Cj > and c§ > such that 

F(s) > c 7 \s\ 2+2q - c 8 , VsGK. (3.14) 

Moreover, (H6) leads to establish further regularity properties for ip, ip t , Uf. This is stated 
in the following 

Corollary 1. Suppose that the assumptions of TheoremUl with (H3) replaced by (H6). 

Then, for every T > there exists a weak solution [u,ip] to Problem P on [0,T] corre- 
sponding to [uq, ipo] such that 

up G L°°(0,T;L 2+29 (ft)), (3.15) 

(ft G L 2 (0, T; V'), if d = 2 or d = 3 and q > 1/2, (3.16) 

u t EL 2 (0,T-X lv ), if d = 2, (3.17) 

and 

'f{tT 3/ \Zlt\>^, (3 - l8) 

2(2v — 3) 

<PteL"(Q,T;V3, if d = 3, 3/2 < p < 2, < q < ^-f (3.19) 

o — p 

where s — ((4 — d)p + 2d)/2p and in (I3.19P the exponent a is given by 

2p(l - f ) 

a - 



(2p-3) -g(3- §)■ 

In two dimensions, as further consequences of (H6), we can prove the energy identity 
and a dissipative estimate, provided that h G £ 2 (0, oo; V£ iv ). Indeed, we have 

Corollary 2. Let d = 2 and suppose that the assumptions of Theorem[J\with (H3) replaced 
by (H6) hold. Then the weak solution [u, ip] to Problem P corresponding to [uq, uo q ] satisfies 

^(«,¥») + i/||Vu|| 2 + ||VAi|r = </»,«>. (3.20) 
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Therefore, (I3.13P with the equal sign holds for every t > 0. Furthermore, if in addition 
h G L 2 (0, oo; V^ iv ), then the following dissipative estimate is satisfied 



where m = (ip , 1) and k, K are two positive constants which are independent of the initial 
data, with K depending on Q, v, J, F, ||/i||L2(- 0)00 .^/. ). 

Remark 7. It follows from Corollary [2] that, in two dimensions, u G C([0, T}; Gdiv) and 
ip G C([0, T]; H). This fact along with the validity of an energy identity suggests that 
the generalized semiflow approach devised in jS] (see also [21]) might be applied to our 
system. If so, one should be able to establish the existence of a global attractor. This is 
one of the issues which will be investigated in a forthcoming paper. 

4 Proof of Theorem d 

The proof will be carried out by means of a Faedo-Galerkin approximation scheme. We 
will assume first that ip G D(B). The existence under the stated assumption on ip 
will be recovered by a density argument by exploiting the form of the potential F as a 
quadratic perturbation of a convex function (see Remark [2]). 

We introduce the family {wj}j>i of the eigenfunctions of the Stokes operator A as 
a Galerkin base in Vdi V and the family {ipj}j>\ of the eigenfunctions of the Neumann 
operator 



as a Galerkin base in V. We define the n— dimensional subspaces W n := (u>i, • • • ,w n ) 
and *5/ n := (-01, • • • , ip n ) and consider the orthogonal projectors on these subspaces in Gdiv 
and H, respectively, i.e., P n := Pw n and P n := Pq, n . We then look for three functions of 
the form 



S{u(t), <p(t)) < £(u , ¥ )e~ kt + F(m)\n\ + K, 



Vt > 0, 



(3.21) 



B = -A + I 



n n n 



U 



n 



k=l k=l k=l 



which solve to the following approximating problem 



((p' n , VO + (Vp(-, <p n ), VV>) = f(u n - + [ (V J * (p n ) ■ VV> (4.1) 
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for every ip G \I/ n and every w G W n , where y?on = -PnV^o and Uo n = P u Uq (primes denote 
derivatives with respect to time). In ( 14. 2 p h n G C°([0, T]; G^) and, on account of (H5), 
we choose the sequence of h n in such a way that h n — > h in L 2 (0, T; V^ iv ). It is easy to see 
that this approximating problem is equivalent to solving a Cauchy problem for a system 
of ordinary differential equations in the 2n unknowns a\ , 6^ . Since i 7 " G C 1 (1R), the 
Cauchy-Lipschitz theorem ensures that there exists T* G (0, +oo] such that this system 
has a unique maximal solution a*™) := (a^, • * ' 5 a n^), := (b\ , ■ ■ ■ ,b^) on [0, T*) 
and aW, b (n) G ^([O, T n *); E n ). 

We now derive some a priori estimates in order to show that T* = +oo for every 
n > 1 and that the sequences of <p n , u n and p n are bounded in suitable functional spaces. 
By using p n as a test function in (14. ip . u n as a test function in ( I4.2p and recalling that 
b(u n ,u n ,u n ) = (see ( 12. ip ). we obtain 

((p' n ,Hn) + (Vp(-,y?„), V/X n ) = / (« n • V/i„)<^„ + / (V J * (f n ) • V/i„ 

z/||Vw„|| 2 = -/(«„• Vp„)y?„ + (h n ,u n ). 
Jn 



1 ^ ii2 

2^' K " 



We now have 



= 11(2 II V^nf + / f(^) - 2&n, J*<Pn. 




J(x - y)((p n (x) - tp n {y)) 2 dxdy + / F(<p n )J. (4.6) 



n Jn Jn 



d /l 

Furthermore, observe that 

{Vp(-,(p n ), VHn) = (-p{-,(p n ),A(i n ) = (-p n , A/i„) = (Vp„, V/i n ), 

where p n := P n p(-,(p n ) = p n + P n (J * <p n ). Summing the first two identities and taking 
the previous relations into account we get 

2~^(ll M «H 2 + \ j f J(x-y)( ( Pn(x)-ip n (y)) 2 dxdy + 2 J F(tprS) 
+u\\ Vu n \\ 2 + || Vp n || 2 + (V(P n (J * tp n )), Wp n ) 

= / {VJ*<p n ) ■ Vp n + {h n ,u n ). (4.7) 
Jn 

Now, it is easy to see that 

||V(P„(J*</v))|| < WB^P^J*^ 

< ||VJ * ip n \\ + || J * ip n \\ < || J||vyi.i||y?„||, (4.8) 
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and that, by means of (H3), we have 

1 
2 




J(x - y)((p n (x) - (p n (y)) dxdy + 2 / F((p n ) 
'nJn Jn 



= \\^(f n \\ 2 + 2 F(<p n ) - (<p n , J * (fin) 

Jn 

> (a + 2c 1 -\\J\\ L i)(p 2 n -2c 2 \n\>a\\(p n \\ 2 -c } (4.9) 
Jn 

where a = 2c\ — \\J\Il 1 > 0. Hence, integrating ( 14. 7p with respect to time between and 
t G (0, T*) and using (" 14.8 j) . ( I4.9j) . we are led to the following differential inequalities 



\u n \\ 2 



aW^nW 2 + / (^||Vm„,|| 2 + ||V/i„|| 2 )(ir < c|| J||^i,i / ||y?„|| 2 dr 
Jo Jo 



+ !l U On|| 2 + 7T / / J ( X ~ y)(<POn{x) - ¥0n{y)) 2 dxdy 



2 



in Jn 



+2 f F(ip 0n ) + ^- [ \\h n \\ 2 v , dr + c 

<M+-*-/ \\hf v , dr + c I y n \\ 2 dr, VtG[0,T n *), (4.10) 
^ Jo dlv Jo 

where c only depends on || </||vpi,i and on \Q\, while M is given by 

M = c(l + ||u || 2 + ||v?o|| 2 + J F(<p )y (4.11) 

Here we have used the fact that, since ipo is supposed to belong to D(B), then we have 
fon ~~ >* V^o m H 2 (Q) and hence also in L°°(Q) (for d = 2,3). Since we have ||tt n (£)|| = 
|a^(i)| and ||<£>n(£)|| — |b (^)l? by means of Gronwall lemma we deduce that T* = +oo, 
for every n > 1, i.e., problem (I4.ip - (I4.5I) has a unique global in time solution, and that 
(I4.10p is satisfied for every t > 0. Furthermore, we obtain the following estimates holding 
for any given < T < +oo 

||«n|U~(0,T;G di „)nL2(0,T;%„) < N, (4.12) 
lkn|U«(0,T;H) < N, (4.13) 
||V// n || L2(0 ,T;i/) < AT, (4.14) 

where 

N = cM 1 / 2 + c\\h\\ LHo , T . yL) , 

with c now depending also on T and on v. From (14.41) . (I4.14p and recalling (ll.Sp we now 
deduce an estimate for ip n in L 2 (0,T; V). We have 

(/i n , -A(p n ) = (V/i„, Vifn) = {-A(f n , wp n + F'((f n ) -J*<p n ) 

= {Vcpn, aV<p n + y?nVa + F"((p n ) Vy? n - V J * (p n ) 
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> c ||Vv?„|| 2 - 2||VJ||ii||V^ n ||||^„|| 

> j\\V<Pn\\ 2 -ky n \\ 2 , (4.15) 
where k = (2/co)||VJ|||i and where we have used (H2). Since 

(V/i n ,V^) < ^||V^ n || 2 + -||V/i n || 2 , 
4 c 



we get 

and dUS}, fOijl . dSll) yield 



||V// n || 2 >|||V^|| 2 -c||^|| 2 , (4.16) 



lbn||L»(0,T;V) < ^. (4.17) 

The next step is to deduce an estimate for the sequence of \x n in L 2 (0, T; V). To this aim 
we first observe that (H4) implies that |-F'(s)| < c|F(s)| +c for every sGi and therefore 
we have 

/ A*n = \{Hn, 1)| = |(a^n + F'fon) ~J*<Pn,l)\ 

< [ \F'(<p n )\+c<c f \F{cp n )\+c<N, (4.18) 
Jn Jn 

since we have \\F(<p n )\\ L oor 0tT . L irQ^ < N due to (14 .7p (integrated in time between and 
t G [0, T]) and (14.91) . We have also used the estimates (14. 12j) - (14. 14[) . Hence, by means of 
the Poincare-Wirtinger inequality, from (14.141) and (14. 18p we get 

\W\\l*(o,t;V) < N. (4.19) 
We also need an estimate for the sequence {p(-,f n )}- From (H4) we immediately get 

||p(-,OIU? < (c\\a\\ L oo\\<f n \\ + \\F'(<p n )\\ LV ) <c((^|F(^ n )|) 1/P + l) <JV, 

and hence we have 

\\p(;<Pn)\\L°°(0,T;L P (n)) < N. (4.20) 

The final estimates we need are for the sequences of time derivatives u' n and (p' n . Let 
us start from the sequence of u' n . Equation (14. 2 p can be written as 

u' n + vAu n + P n B(u n , u n ) = -P n {(p n Vp n ) + P n h n (4.21) 

We now have, for d — 3, by using Sobolev embeddings, interpolation between L p spaces 
and gZESD 

\\P n ((PrS7 p n )\\ V ' Hv < c||v5„,|| L 3||V/i„,|| < c||yJ n || 1/2 ||yJ n ||^ 6 2 ||V)U n || < iV 1/2 ||^„||y 2 ||V/i„,||. 

(4.22) 
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Therefore, thanks to (I4.14p and f)4.17p . we get 

|| Pn{<PnVlin) \\l^(0,T^J < N 2 . (4.23) 

For the case d — 2, by means of Gagliardo-Nirenberg interpolation inequality in dimension 
2 we have, for every < 7 < 1 



\\Pn{^rN ^n)\\v^ iv < c\\ fn || £2+ 7 /(l- 7 ) || V/i n || 

<cii^ir (i - 7)/(2:7) ii^iif 2 - 7 ^/inii 

<iV 2 M/M||^|| 7/M ||V//„||, (4.24) 

so that (gUD and (HT7|I yield 

ll^nC^V^II^^r^j < N 2 . (4.25) 

Moreover, we have ||.Aii n ||y = ||«n|| v div , while the treatment of the term P n B(u n ,u n ) is 
classical and, by means of ( 12. 2 p and f)2.3p we have 



\P n B(u n) u n )\\ v ^ v < c\\u n \\ 1/2 \\u n \\y* v , ford = 3, (4.26) 
\P n B(u n ,u n )\\ v > < c\\u n \\\\u n \\ Vdiv , ford = 2. (4.27) 



Hence, by using (jj^SD , (14251) and flCTj) . (g2ZD, and recalling that P n G ^(V^V^J, 
which implies that 

\\PnK\\l?(0,T;V^ v ) < c(l + ||/i|| £2(0,2^)), 

from (14.211) we obtain 

IK I Il^o.t;^ j < L , for d = 3 > ( 4 - 28 ) 
\W n \\ L 2-, (0i T;V> iv ) <L, V 7 G (0, 1), for d = 2, (4.29) 

where L = N 2 + N. 

In order to derive an estimate for the sequence of tp' n , we aim to take the test function 
ifj G V s in (14.11) . where s > 2 is such that G H S ~ 2 (Q) ^ L p '(fi) (p' is the conjugate 
index to p). Since IP -2 LP* , where p* = 2d/ (d + 4 — 2s), we see that it is enough to 
take 

. > (4 -^ + 2d . (4.30) 

Let us now decompose ^ as 

where = fnV' = ]CLi(V^fc)^fc e ^« and = (/ - P n )^ = Y^^n+ii^^^k e 
(recall that and are orthogonal in all the Hilbert spaces V r , for every < r < s), 
and notice that we have, due to (I4.20p 

l(v P (-,^n),v^)| = !(>(•, o,a^)| 
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Furthermore, it is easy to see that 

/ (v J * <p n ) ■ Wi 

Jn 



< N\\A^\\ LP , < N\\,/>i\\v. 
<NU\\v.. 



< c||VJ|| L l||^r 



< N 



(4.31) 



(4.32) 



As far as the first term in the right hand side of (14. ip (written with ip = ipj) is concerned we 
notice that V^j G /P _1 (fJ). Therefore, when 1 < p < d/(d-l) and s = ((4-d)p+2d)/2p 
orp = d/{d-l) ands > ((A-d)p+2d)/2p = (d+2)/2, due to the embedding H s ~ l m> L°°, 
we have 



/ (U n ■ V^/)v?n 

Jn 



< c||u n ||||^ n | 



(4.33) 



When p = d/(d - 1) and s = ((4 - d)p + 2d)/2p = (d + 2)/2, due to the embedding 
H s ~ l L q for every 1 < q < +oo and interpolation in LP spaces, we have, for every 
r > 2, that 



Jn 



< C\\u r 



V s WuWl^/^- !) 



< c\\Ur 



V s Wn\ 



II 1 1 2/-r- 
VillVnlly • 



Finally, in the case d = 3, when 3/2 < p < 2 and s = ((4 — d)p + 2d)/2p 
due to the embedding iP -1 ^ L 3p// ( 2p ~ 3 ), we obtain 



(4.34) 
(p + 6)/2p, 



< cIIut 



n 

< c||M n |||' 

<N 3 ^\\ v ^ n \\$ p - 3)/p . 



v\\u> Il (3 - P)/P llw ||( 2p - 3 )/ P 



Collecting fl4.3ip - fl4.35p . from (14. ip (written with ^ = ^ij) we then get 

(4 - d)p + 2d 



|^IU°°(0,T;F s ') < L if 1 < p < d', s 



2p 



<Pn\\L°°(0,T;Vl)r\Lr(0,T;V' ) < £ if £> = d', S > 



d + 2 

, r > 2, 

2 ' - ' 



where d' = d/(d — 1), while in the case d = 3, if 3/2 < p < 2 we find 

p + 6 



|< / 9 n || L 2p/(2p-3)( 0iT;V //) < L, S 



2p 



where L = N + iV 2 in all cases. 

From the estimates (jO^ - dPS} . (Q7I . (0511 . (OD]l . ( JOg) , (@22J, © 
and on account of the compact embeddings 



(4.35) 

(4.36) 
(4.37) 

(4.38) 



L 2 (0,T;\/) niZ^O.T;^') L 2 (0,T;#), 



14 



L 2 (0,T; V dlv ) n W^(0,T; V^) ^ L 2 (0,T; G div ), Vg > 1 
we deduce that there exist 

u G L°°(0, T; G dlv ) n L 2 (0, T; V^ w ), (4.39) 

<p G L°°(0,T; if) fl L 2 (0,T; V), (4.40) 

/iGL 2 (0,T;\/), (4.41) 

p G L°°(0, T; L p {Vl)) (4.42) 



with 



^GL 4/3 (0,T;\/jJ, if d = 3, (4.43) 
Ht eL^(0,T;4), V 7 G(0,1), ifd = 2, 



and 



tpteL^T-X), ifl<p<d',s= {4 d) 2 l + 2 \ (4.44) 

r7 _|_ 

^er^TiKjn^o.Tiyy, if P = d', s > , r > 2, 

G L*/(*- 3 )(0,T; if d = 3, 3/2 < p < 2, s = 

2p 

such that, for a not relabeled subsequence, we deduce 

u n — 1 u weakly* in L°° (0, T; , (4.45) 

u n ^u weakly in L 2 (0, T; V^), (4.46) 

u n -^-u strongly in L 2 (0, T; G div ), a.e. in O x (0, T), (4.47) 

u' n ^u t weakly in L 4/3 (0,T;1/^J, d = 3, (4.48) 

u' n - «t weakly in L 2 ^(0,T; l^J, V 7 G (0, 1), d = 2, (4.49) 

V?„ <p weakly* in L°°(0, T; if), (4.50) 

uo n ->> up weakly in L 2 (0, T; 1/), (4.51) 

y? n ->■ <p strongly in L 2 (0, T; H), a.e. in x (0, T), (4.52) 

p n fj, weakly in L 2 (0, T; V), (4.53) 

p(.,<p n )^p weakly* in L°°(0,T;i/(ft)), (4.54) 

and 

<p' n ^<p t weakly* in L°°(0,r;K), (4.55) 
if 1< p < d', with s = ((4 - cf)p + 2d)/2p, 

(ft weakly* in L°°(0, T; V,'), weakly in L r (0, T; V^), (4.56) 

2 
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if p = d', with s > (d + 2)/2 and r > 2, 

tp' n ^<p t weakly in L 2 ^- 3 )(0,T;KO, (4-57) 

if d = 3 and 3/2 < p < 2, with s = (p + 6)/2p. 

We can now pass to the limit in (j4.ip - (|4.5p in order to prove that the functions u and 
p> yield a weak solution to Problem P in the sense of Definition 1, i.e., u, p, p and p 
satisfy fll.Tj) . (I3.9P and ( l3.10p . f )3.1ip . (I3.12p . First of all, from the pointwise convergence 
(I4.52p we have p(-, ip n ) — > a<p + F'((p) almost everywhere in Q x (0, T) and therefore from 
f)4.54p we have p = cup + F'(<p), i.e. (13. 9p . Moreover, since p k = Pk(p(-, <Pk) — J * <Pk), we 
have, for every v G ^ n and every k > n (n is fixed) 

{p k {t),v) X {t)dt= [ {p{;ip k )-J*ip k ,v) X {t)dt, VxGC °°(0,T). 

By passing to the limit as k — > oo in this identity and using the convergences (I4.53p . 
(I4.52p (which implies J * <pk — >■ J * ip strongly in L 2 (0, T; y)) and (I4.54p . on account of 
the density of {^ n } n >i in H we get p = p— J*ip = ap + F'(<p) — J * <p, i.e. ( II. 7p . In 
particular we obtain p G L 2 (0,T; V). 

The argument used to recover (I3.10p and (13.111) by passing to the limit in (14. ip and (14. 2 p 
of the approximate problem and by exploiting the above convergences is standard and we 
only limit ourselves to give a sketch of it. We multiply ( 14. ip by x an d (14. 2 p by oo, where 
X, u> G C^°(0,T) and integrate in time between and T. Due to the above convergences 
we can pass to the limit in these equations. In particular the term (Vp(-, <p n ), V^) can 
be rewritten as (p(-,p n ), —Aip) and (14.541) is used. We also recall that in the nonlinear 
term b(u n , u n , w)oo we exploit the strong convergence (14.471) to pass to the limit. The 
limit equations thus obtained hold for every ip G \l/ n , every w G W n (where n is fixed) and 
every x, oo G C^°(0,T). The density of {^„} n >i and {W„}„>i in V s and Vdi V , respectively, 
allows us to conclude that u, <p, p and p satisfy (13.101) for every ip G V s and (13.111) for 
every v G V^ v . Furthermore, observe that (I3.10P can be written in the form 

(<p t , iP) = -(V/i, ViP) + (it, pViP), (4.58) 

and consider the contribution of the transport term in (14.581) . In the case d = 3, by 
arguing as in (I4.22p we have 

\(u,pViP)\ < TV^IlVulllMl^llVV'll, (4.59) 

while, in the case d = 2, by arguing as in (I4.24p we have 

\(u,pViP)\ < iV r2 ( 1 -' 5 )/M||V«||||^||t/ M ||VV'|| ) (4.60) 

for every 5 G (0, 1). From (14.591) and (I4.60p we deduce that ptif) can be continuously 
extended to V for almost any t > and from these equations and (14.581) we also infer 
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that 

(fit G L 4/,3 (0, T; V'), if d = 3; <p t G L 2 ~ s (0,T;V'), V<5g(0,1), if d=2. 

We hence get (13.61) . (13.71) and furthermore, (I3.10p and (I4.58P hold also for every ip G V. 

Finally, in order to get (I3.12p . it is enough to integrate (I4.ip . (14. 2p between and t 
and pass to the limit for n — > oo by using the weak convergences above. By integrat- 
ing between to and t we prove the weak continuity of u and <p in Gdiv and H, respectively. 

We now prove that the energy inequality (13.131) holds for the weak solution [u, ip] 
corresponding to the initial data no G Gdiv and tpo G D(B). To this aim let us first 
observe that, for almost any t G (0,T) and for a not relabeled subsequence we have 

u n (t) ->■ u(t) strongly in G div , (4.61) 
Pnit) (p{t) strongly in H and a.e. in Q. (4.62) 

and that, by means of (H3) and of Fatou's lemma we have 



F(<p(t)) <liminf / F(<p n (t)). (4.63) 
In addition, it is easy to see that 

P n {J*Vn)^J*V mL 2 (0,T;V), (4.64) 

as a consequence of the convergence J * <p n — )■ J * (p strongly in L 2 (V) and of the fact 
that P n G £(V,V). Hence, by integrating (14.71) between and t, and by passing to 
the limit using (14.611) - (I4.64p . the weak convergences (14 .46 p . (I4.53P and the weak lower 
semicontinuity of the norm, we immediately get (I3.13p . 

In order to complete the proof of the theorem we now assume that uq G Gdiv and that 
fo G H such that F(<p ) G For every k G N let us define ipok G D(B) as 

<p ok := (i + ^-B^j (p . 

Since B is maximal and monotone, we have ip^k <Po in H. Let [uk, <Pk] be a weak 
solution corresponding to u and <y9 0fc and satisfying (I3.2p - (l3.12p . We know that [uk,<Pk] 
satisfies the energy inequality (I3.13P for each k, on the right hand side of which we need 
to control the nonlinear term that, by virtue (|3.1|) . can be written as 

/ F( Vok )= [ G(^ 0fc )-^||^ 0fe || 2 . (4.65) 

To this aim we multipy the equation ip k — ip = —^Bip k by g((pok) in L 2 (Q), where 
g = G'. We obtain 

fi , (v ? ofc)(v ? ofc - (fo) = - j I g{Vok)Bip ok 



17 



= -r f g'{<Pok)\V(pok\ 2 - j; j g(<Pok)<Pak<0, (4.66) 

since g is monotone nondecreasing and we can suppose that g(0) = 0. Therefore, due to 
the convexity of G we can write 



G(ifok) < / G(<p )+ / g(<p ok )(<Pok - <Po) < / G(^ )- (4-67) 
n Jn Jn Jn 

Hence, on account of (I4.65P and (14.671) we get the desired control and from (I3.13p . written 
for each weak solution [u k ,p k ], by means of (H3) and of Gronwall lemma, we deduce 
the estimates (I4.12p . (I4.13P and (I4.14p for u k , (pk and V//jt, respectively. By taking the 
gradient of \i k = a(p k — J * p k + F'(<p k ), multiplying the resulting relation by Vp k in 
L 2 (Q) and using (H2) we recover the control of the gradient of (p k from the gradient of \x k 
(see (14.161) ) and therefore, for ip k we get the estimate (14.171) . Moreover, arguing as in the 
Faedo-Galerkin approximation scheme above we get (I4.19P and (14.201) for \i k and p(-,<p k ), 
and (14.281) . (I4.29P and (I4.36p - ()4.38p for the time derivatives u' k and <p' k , respectively. By 
compactness we hence deduce the existence of four functions u, <p, \i and p satisfying 
( I4.39p -( l4.44p such that the convergences (I4.45p - ()4.54p hold. By passing to the limit in 
the variational formulation for [u k , tp k ] it is immediate to see that [u, <p] is a solution 
corresponding to the initial data Uq and Pq. This completes the proof of the existence of 
a weak solution when uq £ Gdiv arid (po £ H such that F(po) £ L 1 (f2). 

Finally, the energy inequality (13.131) for the solution [u, <p] can be obtained by passing 
to the limit in the energy inequality (13 . 13[) written for each approximating couple [u k , p k ], 
using the weak/strong convergences (I4.45[) - (I4.54[) and Fatou's lemma, in a similar way as 
done above for the Faedo-Galerkin approximate solutions (see fl4.6ip - fl4.63j) ) . In partic- 
ular, on account of (13.11) . when we pass to the limit in the nonlinear term on the right 
hand side we have, by (I4.65p . (14.671) we infer 

limsup / F(p ok ) < / G{p ) ~ tt\\Po\\ 2 = / i^o)- 
fe->oo Jn Jn ^ Jn 

The proof of Theorem [T] is now complete. 

Remark 8. If we compare estimates (I4.28P and (I4.29P for the time derivatives u' n in the 
case d = 3 and d = 2, respectively, with the analogous estimates that hold in the case 
of the local Cahn-Hilliard-Navier- Stokes system (see, e.g., [ID]), we see that in the case 
d = 3 we obtain the same time regularity exponent 4/3 for both the local and nonlocal 
systems. However, in the local system we can estimate p n in L°°(0,T; V) so that, in two 
dimensions we easily get the exponent 2. For the nonlocal system, this possibility seems 
out of reach since we can only estimate p n in L 2 (0,T; V). Also, for the same reason, the 
transport term in the Cahn-Hilliard equation is less regular so that the bound on ip' n is 
weaker in comparison with the analog for the local system. 
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Remark 9. We point out that energy inequality ( 13 . 13[) can be written in an alternative 
form, provided that a suitable condition holds. Indeed, suppose that 



Po = 0, (4.68) 

n 

and that c (see (H2)) and J are such that 

where Cp is the Poincare-Wirtinger constant in the inequality 

IMI < Cp||V^||, \/pE Vs.t. <p = 0. 

Jn 

Then, we can get the following control of the gradient of p by the gradient of [i 

||V/i|| 2 >/3||V^|| 2 , (4.70) 

where ft = (c - 2C P \\ VJ|| L i) 2 (compare (ETTOl with ([US]) ). Indeed, by taking the 
gradient of fi = ap — J * p + F'(p), multiplying the resulting relation by Vp> and using 
(H2) we have 

^l|V^|| 2 + ^=||V/i|| 2 >(V/i,V^) 
>c ||V^|| 2 -2||VJ|| L i|M|||V^|| 

> (c - 2C P ||VJ|| L1 )||V^|| 2 = v^llV^f, (4.71) 

whence f)4.70p . Therefore, as a consequence of (13.131) . for the weak solution [u,ip] of 
Theorem [T] the following energy inequality is satisfied as well 

' f "u{t)\\ 2 + \ [ [ J(x - y)(p(x,t) - p(y,t)) 2 dxdy + 2 [ F(p(t)) 



n Jn Jn 



2 V" v /n 2 
+ l\u\\Wu(r)\\ 2 + P\\Wp(r)\\ 2 )dr 



o 



<\{\\M\ 2j r \ I I J(x-y)(po(x)-p (y)) 2 dxdy + 2 [ F(p ] 
A v A Jn Jn Jn 



+ / (h(r),u(r))dr. (4.72) 



o 



We recall that Cp can be estimated for many important special classes of domains (cf., 
e.g., [28J ) . For example, if Q is convex we can take Cp = diam(fi)/7r and there exist 
convex domains for which this constant is optimal (see [5]). 
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5 Proofs of Corollaries U and [2] 



Proof of Corollary [7J Recalling (I3.14p and repeating the proof of Theorem [TJ in place of 
P~gD we have 




J(x-y)((p n (x)-(p n (y)) dxdy + 2 / F((p r 
>n Jn Jn 



1 
2 

= ||v/a^ n || 2 + 2 I F{ip n ) - (ip n , J * ip n ) 
Jn 

> J {{a - \\J\\l^I + 2c 7 \<p n \ 2+2tl ) - 2c 8 \n\ > c\\<p n \\liX iQ) - c, (5.1) 

and this estimate, by integrating (14. 7p as done above, allows to control the sequence of 
tp n and yields (I3.15p . All the other estimates for ip n , u n , n n and p(-, (p n ) established in the 
proof of Theorem [T] still hold. The only estimates that can be improved are the ones for 
u' n and <f' n . Indeed, for d = 2, in place of (I4.24p we can write 

\\P n (^rN^n)\\v^ v < c||</? n ||i2+2« (n )||V/in|| < iV||V/i n ||, 

and hence we can control the sequence of P n (ip n Vfi n ) in L 2 (0, T; V^ iv ). This control, 
combined with the control for the other terms in (I4.2ip . yields (13.171) . Furthermore, 
as far as the sequence of <p' n is concerned, we can improve estimates fl4.33p -( l4735|) by 
arguing as in the proof of Theorem [T] and by considering the following cases. Choosing 
s = ((4 - d)p + 2d)/2p (cf. fOOj) and flOTjl ). when 1 < p < d' = d/{d - 1), due to the 
embeddings _£P _1 (fi) <— >■ L°°(fi) (if 1 < p < d') or H s ~ l (Vl) L r (f2) for every r < oo (if 
p = d'), we have 

(u n -VtlJl)ipn < c\\u n \\ \\(p n \\L^i(n) \\lp\W e < N\\lp\\v s - ( 5 - 2 ) 

Jn 

The same estimate also holds for the case d = 3 when 3/2 < p < 2 and q > 2(2p — 
3)/(6 — p), where here we use the embedding H 8 ' 1 ^) m> L 3p ^ 2p -^(Q) and the fact that 
6p/(6 - p) < 2 + 2q. Finally, when d = 3, 3/2 < p < 2 and < q < 2(2p - 3)/(6 - p) we 
have 

/ (u n ■ Vif)i)<p n 
Jn 

<? n\U, \\\\m 1 1 2(3— jj) (H-«)/iaC2— 1 1 , _ II (4p-6g+J*?-6)/p(2-g) ,, , 

<iV||^||?^ 69+W - 6)M2 " 9) ||V'lk s - (5.3) 

Hence, on account of (15.21) and (15.31) . from (14.11) (written with if) = if)j) we deduce (13.181) 
and (I3.19P . The improved regularity (I3.16P for tpt can be obtained by estimating the term 
(u, ipVip) in (I4.58P for the case d = 2 as 

\{u,cpViP)\ < C ||V U ||||^|| L2+23(n) ||VV|| < N\\Vu\\\\ViP\\, 



< C||« n ||||^ n || £ ,6p/(6- P )(Q)||V'||v s 
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and for the case d = 3 and q > 1/2 as 

IOwWOI < c||w|| L 2(i +1 /, )(n) ||^|| L2+2s( n)||V^|| < iV||V«||||VV'||. 

□ 

Proof of Corollary [H For d = 2 the regularity properties (13.16P and (13.17P allow us to 
deduce the energy identity for the weak solution. Indeed, in this case we can take and 
v — u(t) in (13. lip and tp — //(r) in (I4.58p . sum the resulting equations and then integrate 
with respect to r between and t. When we consider the duality product ((p t ,fi), we are 
led to the duality (<p t , F'(<p)) which can be rewritten by taking into account that F\ip) = 
g(<p) — cftp, with g G C 1 (IR) monotone increasing. Now, introducing the functional Q : 
H R U {+00} defined as Q(<p) = J Q G(<p) if G(<p) G L l {Vt) and <?(y>) = +00 otherwise, 
we have (see j6j Proposition 2.8, Chap. II]) that Q is convex, lower semicontinous on 
H and £ G dQ(<p) if and only if £ = G'(<p) = g(<p) almost everywhere in Q. In view of 
(I3.16P and of the fact that g(<p) G L 2 (0,T; V), we can use [T51 Proposition 4.2] and get, 
for almost any t G (0, T) 

(<Pt,F'(<p)) = (<pt,g{<p)) - a*(<p t ,ip) = j t (p{v) - yIMI 2 ) = J t J/^)- 
Therefore, on account of this identity, from (13. lip and (I4.58P we obtain 

(|H| 2 + H^ll 2 - (<p, J * (f) + 2 J F(<pj) + u\\Vu\\ 2 + || V/i|| 2 

d_ 
2~dt 



1 d 
2~dl 



' / M| 2 + ^ / / J(x-y)( V (x)-^(y)) 2 dxdy + 2 f F(<pj) 



q Jn 

2 



+i/||Vu|r+ ||V/i|| 2 = (h,u), (5.4) 

Hence we get (I3.20p . Furthermore, by integrating between and t we get the energy 
identity in integral form, i.e, ( I3.13P holds with the equal sign for every t > 0. 

In order to obtain (I3.2ip . let us multiply equation \i = cup — J * <p + F'(ip) by (p in 
L 2 (Vl). We obtain 

{H,V) = \( [ J(x-y)( V (x)-^y)) 2 dxdy + (F'^), V ). (5.5) 
A Jn Jn 

Now, observe that, due to (3.18) and to the convexity of G we have 



and hence 



F(0) > F(s) + -s 2 - (F\s) + a*s)s 



F'(s)s>F(s)-^-s 2 -F(0). 
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Therefore, from (15.51) we get 

1 
2 



<Ji,<p)> \ I I J(x-y)((p(x)-(p(y)Ydxdy^ I F(<p(t)) 




¥>f-c. (5.6) 



2 

Setting ~fi= p| / n A* an d suppose (<p , 1) = for simplicity. Then we have 

(p,<p) = (/i-~p,(p)< C p ||V/x|||M|, 
and then, by means of (H6), from (15.61) we have 

J{z- y ){<p(x)-<p(y)) 3 dxdy+± f i%) + | f W\ 2+2q - c, 




n jo, 



Therefore, we deduce 

Iff J(x-y)( V (x)-<p(y)) 2 dxdy+^ f F(ip) < || V/if + c 10 

and hence 

^K^) < en (^|| Vw|| 2 + || V/if) +c 10 , (5.7) 

where Cn = max(l, 1/2AjJ/), Ai being the lowest eigenvalue of the Stokes operator A. 
We point out that all constants only depend on the parameters of the problem and are 
independent of the initial data. Now, by virtue of (13.201) and (I5.7P we have 

±E{u^) + kE{u^)<l + ^\\h\\l L , (5.8) 
where k = l/2cn and I = cio/cn. By means of Gronwall lemma we hence deduce 

£(u(t), <p(t)) < £(uq, <p )e- kt + K. (5.9) 

with 

If m := (ifo, 1) 7^ 0, observe that if [u, if] is a weak solution with data [uq, ipo] for the 
problem with potential F, then where tp = ip — m is a weak solution with data 

[wo, <Po — m] for the same problem with potential F given by 

F(s) := F(a + m) - F(m). 

By relying on (15. 9p satisfied by the solution [u, ip\, we easily get (I3.2ip . □ 
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